# Call external software to process data (clean, remove rayleigh, scale) 
Drint "Calling external processing software . \n" ; 

system ("process_modis") == 0 or die "Problem with process_modis : $! , 



##SNIPPET FROM process_modis_color . f 90 PERTAINING TO DE** 

#i 1) Atmospheric Correction (snippet) 

#' airmass = 1.0/mu_sun + 1.0/mu_sat 

# trans_sat = exp (-1 • 0*ray_taus (i) /mu_sat) 

# 
# 

I . _ / ,o n/^n + mu sat) + (2 . 0/3 . 0 - mu_sat) *trans_sat) / & 



s^ : ^:j^^f5s%v---- * -pel..*— <« > 

# trayu = ( (2.0/3.0 + mu_sat) + (2 . 0/3 . 0 

# ( 4.0/3.0 + ray taus (i) ) . 

J trayd = ( (2.0/3.0 + mu_sun) + (2 . 0/3 . 0 - mu_sun) *trans_sun) / & 

{ ^•! / ^^ y ^l 1 Uyunr.^«p(-1.0^™..M«»2(i) -wv.i))) ) 
J Proc chandat =-tro/(i.0 + ra y sa(i)*tro) ! units of reflectance [0-1] 
#» 2) Truncate negative numbers to 0 f if they occur 
Si truncate numbers greater than unity to unity j they occur 

# WHERE (proc_chandat < 0.0223) proc_chandat = 0.0223 

# WHERE (proc_chandat > 1.0) proc_chandat - 1.0 
#i 3) Scale the processed data if applicable 

#* print *, "Scaling the data by logJLO" 

# proc_chandat = LOG10 (proc_chandat ) 
#! 4) Dust Enhancement 

jt nrint- * "Creatinq dust enhancement" rt ^ 0 -->* \ 

I Sere (' ( (chanO^ dat + chan02 dat) > 0.0) .AND. (chan02 dat > 0 0223) ) 
t duSJnhance = 7chan02_dat-chan03_dat ) / (chan02_dat + chan03_dat) 

# ELSEWHERE 

# dust_enhance = -0.40 

# ENDWHERE 

#i 5) Truncate between desired values 

# WHERE (dust_enhance < -0.40) dust_enhance - -0.4U 

# WHERE (dust_enhance > 0.15) dust_enhance =0.15 
###END SNIPPET### 

# 

# Dust Enhancement 

# 

if ($have_ir) { 

print "Creating a V'dust enhanced\" image: $ {OUTNAME} . dustlKM. $SECTOR. x. tif f \n 
### ' 

# Water component 

$ $ $ \ if 

print "Doing ocean component ... \n ; 

$min_blue_scl = -1.45; 
$max__blue_scl = 0.0; 
$min_green_scl = -1.45; 
$max_green_scl = 0.0; 
$min_de_scl = -0.40; 
$max_de_scl = 0.15; 

system ("emathp expr_ V ars=- ref_4 69 ref_555 DE' num_ex P rs=3 ". 
"yl_expr='land ? $min_blue_scl : xl ^ ^ 
"y2_expr='land ? $min_green_scl : x2 ' ". 



r ' * 



» Y 3 expr='land ? $min_de_scl : x3' ". n 

Red_ocean-J. units= var _ types= , f loat floa t float' scale_of f sets- scale_f actors= 
" re flec_assem_ $P id tmp_ocean_ $P id") == 0 or die "Problem with emathp: 



$!"; 



system ("imscale include.vars- Blue_ccean. image_colors=250 bg_color=250 
bad.color^^^^^^ under _ color=0 max _width= $ SAMPLES max_height=$LINES 
fiX ed_size=yeI f ». tor=i ^ resample=no est _ range =no min_value= $ min_blue_scl 
max value=$max_blue_scl invert_sca le-no ". dig 

e syrtem m rexp $ im"image format=tiff image_var=Blue_ocean image_colors=250 
^-'"ISrSS^S.wln---- t mP _blue_ $P id blue_ocean.tiff_ $ pid») == 0 or 
diS • "Problem occurred during expim of blue.ocean. tif f_ $P id creation: $ ! «; 

system ("imscale include_vars=. G reen_ocean - image_colors=250 bg_color=250 
bad.color^™.^^^^ under _ color=0 max _ wid th=$ SAMPLES max_height= $ LINES 
fixed_ S ize z =yeI f « ctor=i ^ resample=no est _ range =no min_value= $ min_green_scl 
max value=$max_green_scl invert -jcale-no ". _ Qr 

" ..north_up=yes ^-^VimcSS'S Si2 ocean-component of dust 
"Problem occurred during imscale or green 

enhancement: $!"; ^^..Hff imaae var=Green ocean image_colors=250 

system ("expim image_f ormat-ti£l image_va 
color_palette=black-white "• $ id green ocean. tiff_$pid") == 0 

"draw_indexes= draw_names- tmp_green_9p-L" y 

° r "Problem occurred during expim of green_ocean. tiffjpid creation: $ >»; 

system Ci-cale include_vars= • Red.ocean • image_colors=250 bg_color=250 
bad - COl °Cer'color=249 under_color=0 max_width=$ SAMPLES max_height=$LINES 
fi xed_size=yeI f » ctor=i ^ resample=no est _ range =no min_value=3min_de_scl 

max value=$max_de_scl ^ er ^ a l e -^ de $pid") - 0 or die 

" ^^S^-S&l^cS-S-S- («) ocean-component of dust . 

enhancement: $!"; imaae var=Red ocean image_colors=250 

system ("expim image_f onnat-tili image_vai. 
color palette=black-white -. ocean. tiff $pid") == 0 or die 

- "draw indexes= draw_names= ^P-^-^".:: tif f $pi d"creation: $!"; 
"Problem occurred during expim of red_ocean.tirrj?p 

# Make RGB image ocean. ppm_$ P id") == 0 or die 

system ("tifftopnm red_ocean . tif f_5pia _ 
"Problem with tifftopnm call: $. ; 



/U 



r r»> * 



a n tiff $Did > green ocean . ppm_$pid" ) == 0 or die 
system (-tifftopnm green ocean tiff_$ P1 d g _ 

"Problem with tifftopnm call. ?. , ocean. P pm_$pid") == 0 or die 

system ("tifftopnm blue.ocean^tif^pxd _ 

"Problem with tifftopnm call: 5- , $ id „ } == 0 or die 

system ("ppmtopgm red_ocean.p P m_$pid red. _ 

-Problem with ppmtopgm call. S- , n ocean .pgmjpid" ) == 0 or die 

system ("ppmtopgm green ocean. ppm_J>P - 

"Problem with ppmtopgm call. 5^ , ocean . pgm_$pid" ) == 0 or die 

system ("ppmtopgm blue_ocean ppm.Spid _ 

-Problem with ppmtopgm call . $!» $pid blue ocean. pgm_$pid 

system ("rgbStoppm red ocean ;P£- 9 ?£o££ with rgbStoppm call: |£ 

> tmp ocean. ppm_$pid' ) — 0 or , tiff:t mp ocean one . tiff _S P id > 

product: $1"; 



### 

# Sunglint mask 
### 



## ' ,. nt mask for use with dust enhancement imagery... \n , 

print "Making sunglint mask: tor use 
$min_sunglint =30.0; 

fysS rematnp expr_vars=. sunglint ref_85 3 . num_exprs=l ". 

(X l<( $ mi^ 

:: ° s ° e de^rpK:^ ". 

-ndn-good^r 3ave_exprs=< 1 var na m sca ie_factors= "• 

«var _ units= var_types=' float !"^-°== 0 or die "Problem with emathp: 
"reflec_assem_$pid tmp_glmtl_$pid ) 

$! " ; , delude vars='glintmask' image_colors=250 bg_color=250 

system ("imscale include_vars g 

bad - COl0 r o ° ve :- C olor=252 under_color=2 max_width=$ SAMPLES max_height= $ LI NE S 
fiXed - Si ^oof flctor^l real_resam P le=no est_range=no min_value= $ min_green_scl 
m value= $ max_green_scl invert scalane . d ., _ 0 or die 

"north_up=yes " ^ le Y f gre en ocean-component of dust 

"Problem occurred during imscaie ot y 

enhancement: $!"; 

* r-^t-tiff image var=glintmask image_colors-250 
system ("expim image_f ormat-tif f imag _ 

^^?£^M''to»J*-~'**»' tmp_glint_ $P id 

-^p^ of ^ alint - tiff - $pid creation: $! " ; 

# Merge Sunglint Mask and Ocean image one . ti f f_$pid 



imgiuerge: $!"; 
### 

# Land component 



A-3 



print "Doing land component .. An" ; 

#### , _ -1 45; 

$max_blue_scl - 0.0, 
$min_green_scl - 
$max green scl = O.u, 
$min>d_scl = 1.30. 
$max_red_scl = 2.7U, 

$m in_rmblnd = -0.10, orint out=no ref lec_asseia_$pxd , 

$min fmflnd = -2.0; 
$ maxlfmflnd = 2.0, 

$max_rmb_ir diti 
$min ch26 = 5.0, 

ars _, ref469Eef _555 r e f .e 4 5 r e f „853,od 1S „c h2 e 

system ("emathp ^ T V num exprs=10 -. ^ (l0 ,xl) :badval' 

m odis_oh31 ^ di y^ d ^ T bad(x6)? pow 10,X3 P Q 3) :b adval • 
.. y 2-expr='land ^ bad X • P ;bad L - » . 

»v3 expr='land && « bad xo a irlnd) ? §max_ 

^-*^-iroS^^*?S ^id", - 0 or die -.robX- with 

... ■ _ rm 1 



"Ici-J-' 5 ^ — 

emathp: 91"; fmf noOT wv26_norm ir_norm- 

include vars= y^Y 3 -" 
system ("copyvar^ncluQ _ 

overwrite_vars=n - 



$1 !f . 



teid-i - o or die " Problem wl 

, i ■ f f $pid ref lec_assem_$p 1 d ) 

A I II . 



* •„ tooo = 3000.0; fmf n orm wv26_norm 

JSSS scale - 3000.0' re£ 555 y3jS_*« - 

sy s?eS l-er^U &vatton' n« «prs-3 • ,. . 

-~3^g H s Ell i 

.•use_degrees f - exprs ='l 2 3 vai_ fac tars= 
.. min _good= sav _ f scale __ of f se ts= scaled 

Red - land :. v 'a; units= var,types- float £ o. ^ ^ ^ emathp: 

., r eflec_asse m _ $ pid tmp_land_ $P id«) - 

or vars=' R ed land' nu» eKprs-1 
astern ("emathp expr va rs _ E ed_scl' - 

fixed size=yes , resa mple=no est_rang 

fixe - .-zoom factor-1 reai -' sca ie=no "• <n iA") — 0 OE die 

■ problem occurred ^ ima , a _oclors-250 

enhancement: fo rmat-tiff imaae_var . _ „ 

system f ^^hite ". mD landb lue_S P id blue_land. taff _«> 

«lc^»^exe S - <*»-—- ^ land . tiIf _ SP ld creation. 

« - .Problem occurred durln, «- ^ ^ bg _ color=25 0 

sy3 tem (-1— * include,^- Creen_ ^nei^-SUHES 

--^r.color^ nnder_colot=0 _ ^ va _ 5min . g ree».el 

^o^curred durln, -cale , lmage _oolora. 25 0 
oSiette-blaC-blte . 



"draw_indexes= draw_names= tmp_landgreen_$pid green_land.tiff_$pid") == 

° ° r die .. Prob iem occurred during expim of green_land. tif f_$pid creation: $!"; 

system ("imscale include_vars= ' Red_land' image_colors=250 bg_color=250 

bad - COl °«over_color=249 under_color=0 max_width=$ SAMPLES max_height=$LINES 
■fixsd size = ves 11 • 

~ »zoom_factor=l real_resample=nO n est_range=no min_value=$min_red_scl 
max value=$max red scl invert scale=no " . 

max_value^max_ tmp _ lan d_ $pi d tmp_landred_$pid" ) == 0 or die 

"Problem occurred during imscale of red land-component of dust 

"ySerrexpim'image format-tiff image_var=Red_land image_colors=250 

COl ° r - pa i^:1i£;:i i d?aw:names= tmp_landred_ $P id red.land. tif f_ $P id") == 0 or 

^ -Problem occurred during expim of red_land. tif f_$pid creation: $!"; 

### 

# Make RGB image 

"ystem ("tifftopnm red_land. tif f_$pid > red_land.ppm_$pid") == 0 or die 
"Problem with tifftopnm call: $!"; n H ; p 

system ("tifftopnm green_land.tif f_$pid > green_land.ppm_$pid ) — 0 or die 

"Problem with tifftopnm call:. $!"; ^h-i — n or die 

system ("tifftopnm blue_land.tif f_$pid > blue_land.ppm_$pid ) — 0 or die 

l^r,^^-^";^ > red.land.pg^pid", - 0 o, die -Problem 

land.ppmjpid > g ree»_land.p*».Spid« ) - 0 or die 

."pSto^s^aid.Ppipid > bl „ e _ la „d. Pgm _ 5 pi,») - 0 o r die 

"SSr.^S^^.^iO gI ee„_land. P5 n,. 5 pid Mue_land.pgM.Spid > 

tmp land.ppm_$ P id") ==0 

or die "Problem with rgb3toppm call: $! , a^nm — n 

svstem ("convert ppm:tmp land.ppmjpid tif f :tmp_land. tif f_$pid ) — 0 
system ^co^ .^J*^ *- curred during convert of m odis dust enhancement 

product: $!"; 

# Merge Ocean and Land RGB imagery, and metafile image 
system ("imgmerge background_file=tmp_ocean.tiff_$pid 

with imgmerge: $!"; 

system ("legend line_l=\ "$ { satname } -MODIS $mm/$dd/$yyyy $(hhmm)-Z 1km Dust NRL 
MontereyX" k=possible Dust , B lack/White=$ {min_sunglint } -Degree 

Sunglint Zone"" iineV text_height= $ TEXT_HEIGHT text_type=bold line_of f set=0 
sample.offset^V^^ center _ text=n0 outline=yes so lid_bg=yes $ MASTER SLEGEND") 
== 0 

or die "Problem with legend: $!"; 
system ("xinject bg_var_name= -master ■ colors_ P er=2 meta_colors- 5 0 3 0 5 4 



H 



"line_widths='2 1 2' line_types= marker_types= marker_sizes= 
"image_colors=4 $ MASTER $GRID $COAST $ LEGEND meta . tdf_$pid" ) == 0 
or die "Problem with xinject: $!"; 

system ("expim image_format=tif f image_var=master image_colors=4 
"colorj)alette=black-white draw__indexes=\"3 4 5 6\" ". 

"draw__names=\ "violet white blue salmon\" meta.tdf $pid meta. tiff $pid") 
==0 ~ 
or die "Problem with expim: $!"; 

print "Merging composite with meta image. . An"; 

system ("imgmerge background^ ile=$ {OUTNAME} . dust_tmp.tif f_$pid 
foreground_file=meta. tif f_$pid clear_color=black ". 

"output__format=tiff $ {OUTNAME} . dustlKM. $SECTOR. x . tif f_$pid" ) == 0 or di 
"Problem with imgmerge: $!"; . 

# Adjust image quality 
system ("convert -quality 80 tif f: $ {OUTNAME} . dustlKM. $ {SECTOR} . x. tif f_$pid 
"jpeg :$ {OUTNAME} .dustlKM.${ SECTOR} . x. jpg_$pid" ) ==0 
or die "Problem occurred during convert of modis dust enhancement 

product : $ ! " ; 




